^ vr ■' ■ 

/ >_. / " / / j t 

3 


' (' 


N89 - 14 59 6 / 

■-tv*- 

V 
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ABSTRACT. The largest amplitudes of the observed Quasi Biennial Oscillation (QBO) in column ozone are 
found in high latitudes and this must be taken into account in any explanation of the increased deple ion o ozone 
in the southern polar spring during the 1980’s. A QBO in zonal wind, temperature and column ozone has been 
successfully modelled in a two dimensional dynamical/chemical model by the introduction of a parametma ion 
scheme to model the transfer of momentum to the zonal flow associated with the damping ofv " tl <^ ly pr ° paga ln & 
Kelvin and Rossby-Gravity waves. The largest anomalies in column ozone of approximately 20 DU are present at 
ygh latitudes. The equatorial ozone QBO is out of phase with the mid- and high latitude ozone QBO, in good 

agreement with observations. 

1 Introduction , . 

The presence of a quasi-biennial oscillation (QBO) in zonal winds temperature arid column 
ozone in the equatorial stratosphere is well known (Reed 1960, Veryand and Ebdon 1961 V . allace 
r^CllTd Korshover 1978, Coy 1979, Tolson 1981, Hasebe 1983, Naujokat ,1986). The 
largest amplitudes of the ozone QBO, however, are found in high latitudes and this must be 
allowed for in any explanation of the increasing depletion of ozone in the southern polar spring 
during the 1980’s (Bojkov 1986, Garcia and Solomon 1987). We describe a simulation o le 
QBoln a two-dimensional model which extends from pole to pole and includes a 
nhotochemical scheme. We are therefore able to investigate the effects of he equator^ QBO 
on the extra-tropical dynamics of the model and the relationship between the dynamical QBO 

and the distribution of ozone. 

o TVig Model 

The two-dimensional model of Harwood and Pyle (1975,1977,1980) 

The model configuration is essentially the same as that employed by Gray and Pyle (1987a) 
hereafter referred to as GP, except that the improved ‘wide-band’ radiation scheme deve ope y 
Haigh (1984) has been used. This has resulted in a much deeper and stronger Hadley circulation 
in the equatorial lower stratosphere. (Note that, in this respect the present mo ^ el ^ eTS ^ 
that described by Gray and Pyle 1987b). The parametnzation scheme to model the S 
Annual Oscillation (SAO) developed in GP was retained in the present runs but the gravi y 

wave breaking by deluding the momentum deposition associated 

wit^ ther^dly^dissipating ^Kewfn and ^Gravity wave, 

to derive an expression for the mean flow acceleration: Tt - Aexp{ H )R{z)exp{ ( )) 
n( . _ gjfjN and P(z) = f z R(z')dz' (Lindzen and Holton 1968). A is the vertical momentum 

si! at S km) 0 ( 2 ) is the thermal damping rate, N is the Brunt- Vaisala frequency and 
the^onal waveLrier. The equation was s solved with paramet er values ap^oPna^the 

with the faster phase speed c=6C ) ms- occurs a. Uggg* ^‘tv^htlin^Wed to 
gives rise to the 30 km the damping rate 0 ( 2 ) 

wi™i^ob!' the 'slow' damping rate of Dnnkerton (1979) which peaks at approodmately 

9xUT 6 s" 1 at 50 km- below this level a constant value of 0.35x10 s was specified. Although 
2x10 s a ) ii itiv employed in the QBO parametnzation scheme, the model 

irmlude^^vertic^^l^ion^perating^ontne model win?8eld^ with a constant value of = 1.0 

24) 



m 2 s 1 at all latitudes and heights. A gaussian distribution about the equator was applied to the 
forcing associated with the two types of waves, with an e- folding width Yl = (j^)* for the Kelvin 
waves and Yl = ( — ^ )a for the Rossby- Gravity waves where ^ = c (Holton 1975). All other 

symbols have their usual meaning and an overbar denotes a zonal average. The parametrized 
forcing associated with the Kelvin and the Rossy-Gravity waves were both present at all times 
of the modelled year. 
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3 The QBO in Zonal Wind 

The time-height section of the zonal wind at the equator from a model run that included both the 
QBO and SAO parametrizations is shown in figure 1 . Four periods of the QBO are evident. The 
semi-annual oscillation is dominant above 35 km and the QBO is the dominant signal between 
10 and 30 km, in good agreement with observational data. Both phases of the QBO exhibit a 
gradual descent with time (at an average rate of just over 1 km per month); the amplitude of 
the modelled QBO is a maximum at approximately 25 km where the winds vary between 20 
ms" 1 and -18 ms" 1 . Thermal wind balance is maintained against thermal dissipation in the 
model by the development of a meridional circulation with adiabatic heating (cooling) in the 
cold (warm) region associated with upward (downward) motion (Reed 1964, Plumb and Bell 
1982). Therefore, the direction of the induced circulation depends upon the phase of the QBO. 
During a westerly (easterly) phase descending (ascending) motion is present at the equator, just 
below the level of maximum wind shear, with rising (sinking) motion at mid-latitudes. The 
direct forcing of the zonal wund in the model due to the damping of vertically propagating 
equatorial waves is restricted to latitudes less than about 20 degrees from the equator; however, 
changes to the modelled atmosphere, particularly in the amount of the ozone column, occur 
further poleward of this as a result of the induced meridional circulation. 


4 The QBO in Ozone 

Observations of large ozone reductions in high southern latitudes during Spring (Farman et 
al., 1985) make the understanding of the latitudinal variation of the ozone QBO particularly 
important. Garcia and Solomon (1987), for example, have speculated that the QBO is relevant 
to the understanding of the temporal variation of the springtime column ozone in southern polar 
latitudes during the 1980’s. 

The photochemical time constant of ozone deer i ascs with increasing altitude and generally, 
at any given height, with decreasing latitude. With a lifetime of hours in the upper stratosphere 
ozone is expected to be close to its steady state value On the other hand, the ozone concentration 
in the polar lower stratosphere can be far from its equilibrium value since the photochemical 
time scale is long in that region and dynamics exurts a major influence. Between these two 
regions is a portion of the atmosphere w r here dynamical and photochemical time scales are 
roughly comparable. It is the morphology of this transition region that defines the ozone budget 
and distribution. The strength of the lower stratosphere circulation is a crucial parameter in a 
successful model simulation of the ozone distribution. 

Figure 2 shows the time-series of ozone anomaly from the model run (the monthly mean 
has been subtracted from each point). A prominent QBO signal is present at all latitudes. The 
largest anomaly is in high latitudes, reaching 20 Dobson Units in some cases. A phase reversal 
is present at approximately 20 degrees latitude so that high latitude anomalies are out of phase 
with the equatorial anomalies. There is no evidence for in-situ photochemical control of the high 
latitude ozone anomaly. This is not suprising given that the photochemical constant of ozone is 
long there. None of the radicals involved in ozone photochemistry, including OH, H0 2 , NO and 
N0 2 , show significant biennial variations in high latitudes. 

241 







The overall behaviour of the modelled ozone QBO compares favourably with the zonally 
averaged observations of the QBO derived by Hasebe (1983) using SBUV data. In agreement 
with observations the positive ozone anomalies occur at approximately the same time as the 
westerly phase of the equatorial zonal wind QBO at 50 mb (see figure 2). 

The modelled QBO in column ozone arises as a result of the interaction of dynamics and 
photochemistry in the following sequence of events. The equatorial zonal wind QBO give rise 
to a rising and sinking motion at the equator as already described. Hence, during a westerly 
phase the strength of the equatorial ascent of the Hadley cell is reduced and during an easterly 
phase it is enhanced. During an easterly phase the enhanced upward motion lowers the ozone 
column by reducing the lower stratosphere mixing ratios, since the mixing ratio profile increases 
with altitude; conversely, during a westerly phase downward motion increases the ozone column 
by advecting higher mixing ratios into a region where the photochemical lifetime is long. Thus 
during a westerly (easterly) phase a positive (negative) equatorial ozone anomaly will result. 

The induced circulation is also responsible for the high latitude QBO in ozone. During a 
westerly phase of the QBO the induced downward motion at the equator implies an upward 
motion in the subtropics which opposes the Hadley circulation. The latitudinal extent of the 
induced circulation is restricted to the subtropics, however (Plumb 1982), so that it cannot 
directly transfer an ozone anomaly to high latitudes. Nevertheless, the descending arm of the 
Hadley cell extends to mid-latitudes and crosses the transition region from photochemical to 
dynamical control of the ozone distribution. Hence, for example, a suppression of the strength 
of the Hadley circulation during a westerly phase of the QBO results in a negative anomaly 
in mid-latitudes (not just in the subtropics). The ozone anomaly is then transferred further 
poleward by eddy motions (this is opposed by the mean motion and hence the distribution of 
ozone at high latitudes is determined by the balance between eddy and mean motion - see for 
example, Harwood and Pyle (1977) and Haigh (1984)). 

The high latitude column ozone anomaly maxima in the two hemispheres of the run are 
displaced by several months; they tend to occur in spring/early summer and coincide with the 
column ozone maxima. This compares well with observations and arises in the model because 
the phase of the modelled QBO in column ozone tends to reverse during April/May of each year. 
Because of the direction of the mean circulation in the following months the ozone anomaly is 
immediately evident in the southern hemisphere mid-latitudes and this is then transferred to 
higher latitudes via eddy transfer (which is strongest in winter months). The corresponding 
anomaly maximum does not occur in the northern hemisphere until six months later when the 
direction of the mean circulation has reversed in direction. 
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